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Abstract. - The elastic interactions between objects embedded in a nematic liquid crystal are 
usually caused by the average distorsion — rather than by the fluctuations — of the nematic ori- 
entational field. We argue that for sufficiently small particles, the nematic-mediated interaction 
originates purely from the ffuctuations of the nematic director. This Casimir interaction decays 
as , d being the distance between the particles, and it dominates van der Waals interactions 
close to the isotropic-to-nematic transition. Considering the nematic as a polymer solvent, 
we show that the onset of this Casimir interaction at the isotropic-to-nematic transition can 
discontinuously induce the collapse of a ffexible polymer chain from the swollen state to the 
globular state, without crossing the 0-point. 



In 1948, Casimir originally showed that the interplay between electromagnetic quantum 
fluctuations and boundary conditions generates long-range interactions between uncharged 
conducting walls The corresponding energy density per unit area falls off as where 
d is the distance between the two walls. Similar effects can also arise in liquid crystals (H ^. 
Nematic liquid crystals are fluid mesophases composed of rodlike molecules exhibiting a bro- 
ken rotational symmetry along a nonpolar direction, which is characterized by a unit vector 
n called the nematic director Although a nematic material placed between two walls 
imposing identical orientational boundary conditions does not store any elastic energy, and 
hence does not give rise to any mean-field force between the walls, it was shown that the 
interplay between the director's thermal fluctuations and the boundary conditions induces a 
"Casimir" long-range interaction falling off as |3|. Note that Casimir-like interactions have 
been reviewed recently in 

Such Casimir interactions should exist without any distortion of the average director orienta- 
tion, which makes them so intriguing. However, they have not yet been experimentally brought 
to the fore. Indeed, in most studied systems, e.g. macroscopic colloidal particles immersed in 
nematic liquid crystals the interactions mediated by the nematic phase are dominated 

Typeset using EURO-M^ 



2 



EUROPHYSICS LETTERS 



by the energy of the mean-field elastic distorsions of the director field, which decays as 
or depending on the symmetries of the inclusions (|, ||, [oj. We argue here that inclusions 
of size smaller than the molecules's orientational correlation length ^ should experience pure 
Casimir interactions. The latter are found to decay as d~^, and should dominate van der 
Waals interactions in the vicinity of the nematic-to-isotropic transition. Then, considering a 
single flexible polymer chain in a nematic solvent ]lO| , we show that this Casimir interaction 
can induce the collapse of the chain at the isotropic-to-nematic transition. 

In a nematic liquid crystal, the constituent molecules are locally oriented on the average 
parallel to some direction n. Close to the isotropic transition, although the macroscopic degree 
of order can be very weak, the molecular orientations are correlated over "swarms" of size ^, 
which can be as large as 10 to 20 times the molecular dimensions ||^. On length-scales larger 
than ^, these essentially uncorrelated swarms fluctuate about their common orientation n, the 
fluctuations being all the more important than the degree of order is small. To construct the 
nematic order-parameter Q — a traceless tensor describing the amount of uniaxial anisotropy 
of the nematic phase ||^ — one can average over a small volume the microscopic quantity 
q = I (mm — il) made from the molecular orientation m and the identity tensor I. However, 
it follows from the previous discussion that the averaging volume must be of dimension larger 
than ^ for this average to be representive of the macroscopic degree of order. One actually 
defines the director n and the amount of order S by writting the average of q in the form 
Q = |5(nn — ^I). Hence, both S and n are only defined on length-scales larger than the 



Let us now consider a microscopic (neutral) impurity of size less than ^, embedded in 
the nematic phase. On the range of molecular forces (a few molecular lengths), it will 
directly affect the orientation of the nematic molecules, either reducing or enhancing their 
parallelism (see Fig. |l|). Due to the orientational correlations of the nematic molecules, this 
direct influence will spread over ^. However, the nematic director — which is only defined at 
length-scales larger than ^ — will not be affected by the presence of the particle. This is because 
the orientational correlations essentially vanish beyond ^ Such a microscopic impurity 

cannot therefore induce any mean-field distorsion of n, as a macroscopic colloidal particle 
would (the difference being that the latter is physically correlated on a macroscopic scale). 
Nevertheless, a microscopic impurity should affect the local degree of order S, because the 
swarms around it are strongly disoriented. Hence, within the macroscopic Landau-de Gennes 
description of nematics — which is naturally coarse-grained over ^ — molecular impurities are 
simply represented by point-like variations of the degree of order S. Moreover, since the Frank 
elastic constants associated with the distorsions of the director field are proportional to 
such impurities should be modelled as point-like variations of the Frank elastic constants. 

To simplify, we describe the nematic elasticity in the one constant approximation and 
we assume that the director undergoes only small deviations with respect to e^. Setting 
n ~ (rix, riy,!) and considering a spatially varying elastic constant 



which models two molecular impurities located at ri and r2, the system's free energy takes 
the form 



Throughout the paper, all the energies are expressed in units of /cbT. The minimum of (|^) 
corresponds to the uniform state = Uy = 0. Therefore the inclusions do not experience any 
mean-field interaction. The total free energy of the system, however, which takes into account 



cutoff ^. 



K{r) ^ K + Sti 5{r ~ ri) + St2 5{r - r2). 
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Fig. 1. - Snapshot of a nematic liquid crystal in the vicinity of a disordering impurity (black sphere). 
The nematic molecules under the direct influence of the impurity are drawn in dark gray. The 
correlation "swarms", of size 5, indirectly influenced by the impurity are drawn in light gray. Note 
that the orientational fluctuations within the swarm have been diminished for the sake of clarity. The 
director n, which is only defined by spatial averages over distances larger than ^ is not influenced by 
the impurity. In particular, the above +1/2 and —1/2 disclinations of the molecular orientation are 
irrelevant, since their separation compares with ^. However, after coarse-graining on a scale larger 
than ^, the impurity is seen to reduce the nematic degree of order S, thereby inducing a point-like 
reduction of the Frank elastic constants. 



the fluctuations of the director field, depends in a non-trivial way on r = |r2 — ri|. Because 
Tlx and Hy are decoupled and play a symmetric role, the partition function Z, and hence the 
total free energy, can be expressed as a functional integration over n^. only. By performing a 
Hubbard-Stratonovich transformation [p^ , we obtain 

Zi= f PK] exp { - fdr [iif (Vn,)2 + ^ J] 5{r - r,) (Vn^f] } 



hi hi 



(ihidhze 2*1 2 *2 exp | - /"dr [^^(Vn^)^ + i ^ hj • Vn^, (5(r - r^)] |. (3) 



i=i,2 



Then, calculating all the Gaussians integrals yields 

(1 

Z{r) cx det"^ 



I -VVG(r) ^-VVG(O)) 



VVG(O) -VVG(r) 



(4) 



in which the above 6x6 matrix, made of four 3x3 blocks, involves the identity tensor I, and 
the tensor VVG(r). Here, G(r) is the Green function of the operator —KV^. The total free 
energy of the system, consisting of the nematic liquid and of the two impurities, can then be 
calculated as F{y) = — InZ(r). This total free energy can be decomposed as the sum of the 
free energies of individual impurities in the nematic solvent (which we call the ground-state 
free energy), and of the interaction free energy between the two impurities. In the following, we 
consider only the latter, which we denote by Fc(y) and which goes to zero when the separation 
of the particles goes to infinity. The ground-state free energy appears as a constant that we 
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drop. 

To calculate the interaction free energy Fc{r), one has to introduce a short-distance cutoff 
A~^. Note that this cutoff is arbitrary jl^, however it cannot be smaller than the correlation 
length of the nematic solvant, since the free-energy is not defined below shorter distances. 
After regularization involving the cutoff A in the reciprocal space, we obtain (in units of /cbT) 

with SKi{A) = (27r)^A~^(5fci(A). The perturbations Ski{A) depend on the cutoff; this depen- 
dence may be determined from the property that the interaction free-energy, however, should 
not depend on the cutoff. Anyhow, when A~^ is chosen to be of order ^, the perturbation 
of the nematic order due to the impurity is expected to be of the same order as the nematic 
order itself. Thus, from now on, we set A = 27r/^ and we assume \6Ki\ < K. 

We find that the Casimir interaction is attractive when the two inclusions are of the same 
nature, i.e. when SKi SK2 > 0, and repulsive otherwise. Let us compare it with the van der 
Waals interaction, which decays with the same power law. With F^dwir) = —A/Trp^r^, in 
which A is the Hamaker constant and p the inverse of the molecular volume of the impuri- 
ties psi, we obtain 



Fc 27 f 5K ^"^ 



Fvdw 8ttA\3K + SK 



{pey^^^ipey. (6) 



With typically A ~ f (in units of /sbT) and \6K\ < K, we expect p in the range 10 ^-0.25 
(the latter value corresponding to 6K ~ ^K)- Deep in the nematic phase, where p^^ — 10, 
the Casimir and van der Waals interactions should therefore be comparable. However, in 
the vicinity of the nematic-to-isotropic transition, where ~ 10'^, the Casimir interaction 
should be the leading one. There can be dramatic consequences to this fact, as we argue now 
by considering a flexible polymer embedded in a liquid crystal solvent . 

For flexible polymers, the persistence length corresponds to one monomer length and is 
about SA, which is smaller than the correlation length ^ of the nematic liquid crystal. Thus, 
as discussed above, the polymer cannot induce director gradients on the scale of one persistence 
length (as a particle larger than ^ would, e.g. by means of anchoring interactions). Beyond the 
chain persistence length, the monomers orientations are uncorrelated, or, if they are correlated 
[ p^ , this correlation is induced by the surrounding nematic field, and cannot therefore induce 
any director gradient. We conclude that the presence of a non-collapsed flexible polymer in 
a nematic liquid does not give rise to Frank elasticity. In the following, we shall neglect the 
possible large scale anisotropy of the Gaussian conformation of the polymer induced by the 
nematic solvent. 

We suppose that the liquid crystal in its isotropic phase is a good solvent for the considered 
flexible polymer. This property is characterised by a positive excluded volume parameter 
vq, which takes into account the short-range microscopic monomer-monomer, nematogen- 
monomer, and nematogen-nematogen interactions. The Flory theory predicts a collapse of the 
polymer chains when the excluded volume is negative. In the particular case vq = 0, called the 
8-point [|l9), the conformation of the chain is essentially Gaussian. For positive values of the 
excluded volume, chains are swollen as compared to the Gaussian random walk, and for chains 
long enough, the end-to-end distance of the polymer is given, in the Flory picture, by 



, 3 



{R^)2 ^b^v^N^, (7) 
where b is the monomer length and N is the degree of polymerization of the chain. 
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The issue we address now is the effect of the isotropic-to-nematic transition on the chain 
conformation, in the case of a flexible chain with persistence length b < ^. In the ncmatic 
phase, the effective interaction between monomers results now from a bare excluded volume Vb 
that takes into account the short-range interactions, plus the long-range Casimir interaction 
calculated previously. Since at length-scales smaller than ^, the nematic and the isotropic 
phases are essentially indistinguishable from the point of view of the local molecular organ- 
isation (as argued above), we do not expect the isotropic-to-nematic transition to shift the 
value of the bare excluded volume: thus we take Vb = vq. The Casimir interaction, however, is 
defined only on length-scales larger than i^. Hence, to compare the bare excluded volume and 
the Casimir interaction, one must coarse grain the polymer in order that the monomer length 
becomes ^ . This coarse-graining of the chain imposes to rescale the bare excluded volume 
vq. Thus, we redefine the monomer according to 

N N/\ = N', 

b^by = b', (8) 

while the end-to-end distance, which is unchanged, is now given by 

^5'f„;*7V't, (9) 

where is the rescaled excluded volume. Thus, one has 

v'a = vo\\ (10) 

with 7 = 9/5 if we take for v the Flory exponent v — 3/5, or 7 = 2 if we assume that the chain 
has a Gaussian structure at small scale, which corresponds to = 1/2. Note that rederiving 
the end-to-end distance by the Flory argument — i.e. calculating a Gaussian entropy for the 
chain, neglecting correlations to calculate the excluded volume interaction, and balancing 
both contributions — yields the result 7 = 2. Considering now monomers of size 6' = ^, we 
take A = {^/b)~ , and thus the rescaled excluded volume is 

'^o = (J^ " , (11) 

with j/h' = 3 when assuming a swollen state at small scale, or 7/1/ — 4 when assuming a 
Gaussian conformation at small scale. If one considers chains not far from the 0-point, since ^ 
is not order of magnitudes larger than &, we can assume a Gaussian behaviour at scales smaller 
than ^ A rescaling exponent 4 is thus more reasonable than 3, and will be retained for 
the discussion below. 

The chain can now be seen as an ensemble of monomers of size ^, submitted to the bare 
rescaled exluded volume interaction Uq, and to a long-range Casimir interaction characterized 
in our model by SK ~ K. Note that the Casimir effect is pairwise additive at lowest order. 
The effective excluded volume interaction is UofF ~ v'q + vq, with 

., = £*fcM = -|;(5^)V.-o«». (12) 

This Casimir contribution is indeed attractive, and its effect is to reduce the quality of the 
solvent. Adding the two contributions, we obtain 
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Ves^voi^] -af, (13) 
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with a ^ 0.1 for SK ~ K. The collapse condition Vcs < can be written as 

Assuming b/^ ~ 0.1, we find that an abrupt collapse may be induced by the isotropic-to- 
nematic transition, provided the isotropic phase is not too far from the 0-point. For instance, 
with vo/b^ ~ 10~^, chains with degree of polymerization N > 10"*, which are thus swollen in 
the isotropic phase should effectively collapse at the transition. Note that the transition 
from the swollen to the collapse state occurs without crossing the 0-point: this is a consequence 
of the first-order nature of the isotropic-to-nematic transition, which induces a discontinuous 
change of the interactions between the monomers. One could not expect this effect with a 
continuous second-order phase transition. 

In conclusion, we have shown that small impurities in nematic solvents undergo long-range 
Casimir interactions decaying as d~^. For such inclusions, elastic interactions are absent, while 
Casimir interactions dominate over the van der Waals interactions. These Casimir forces are 
attractive for identical particles and can lead to the collapse of flexible polymer chains. They 
could thus be measured by determining the phase diagram of various flexible polymers close 
to the isotropic-nematic transition. 



We thank A. Ajdari and L. Peliti for useful discussions. 



REFERENCES 

[*] New address: Universite de Paris-Sud, Laboratoire de Physique des Solides, bat. 510. 91405 Orsay 
Cedex, France. 

[1] Casimir H. B. C, Proc. Kon. Ned. Akad. Wet. B, 51 (1948) 793. 

[2] Ajdari A., Peliti L. and Prost ,L, Phys. Rev. Lett, 66 (1991) 1481. 

[3] Ajdari A., Duplantier B., Hone D., Peliti L. and Prost J., J. Phys. 11 France, 2 (1992) 
487. 

[4] Kardar M. and Golestanian R., Rev. Mod. Phys. (colloq.), 71 (1999) 1233. 
[5] De Gennes P.-G. and Prost J., The physics of liquid crystals, (Oxford science publications) 
1993. 

[6] PouLiN P., Stark H., Lubensky T. C. and Weitz D. A., Science, 275 (1997) 1770. 

[7] Poulin p., Frances N. and Mondain-Monval O., Phys. Rev. E, 59 (1999) 4384. 

[8] RUHWANDL R. W. AND Terentjev E. M., Phys. Rev. E, 55 (1997) 711. 

[9] Lev B. I. and Tomchuk P. M., Phys. Rev. E, 59 (1999) 591. 
[10] For a review on liquid-crystalline polymers in solution, see Cotton J. P. AND Hardoin F., 

Prog. Polym. Sci., 22 (1997) 795. 
[11] Only the fluctuations associated with the Goldstone mode of n are long-ranged. 
[12] Frank F. C, Disc. Faraday Soc, 25 (1958) 19. 

[13] Chaikin P.m. and Lubensky T.C., Principles of condensed matter physics, (Cambridge uni- 
versity press) 1995 
[14] The Frank elastic constants depend however weakly on this cutoH'. 

[15] Israelachvili J. N., Intermolecular and Surface Forces, (Academic Press, London) 1992. 
[16] De Gennes P.-G., Polymer Liquid Cristal, (Academic N. Y.) 1982. 
[17] Brochard-Wyart F., C. R. Acad. Sc. Pans, 289 (1979) 229. 
[18] Triantafillou M. and Kamien R.D., Phys. Rev. E, 59 (1999) 5621. 
[19] De Gennes P.-G., Scaling concepts in polymer physics, (Cornell university press) 1986 
[20] Constituents of the chain separated by a distance smaller than ^ can be considered as contributing 
to the "charge" &K associated to the corresponding swarm of size ^. 



